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ORIGINAL ARTICLE 

Blood transcriptomic biomarkers in adult primary care 
patients with major depressive disorder undergoing cognitive 
behavioral therapy 

EE Redei^'^'^ BM Andrus^ MJ Kwasny^'^ J Seok'^'^ X Cai^'^ J Ho^'"^ and DC Mohr^'^'"^'^ 

An objective, laboratory-based diagnostic tool could increase the diagnostic accuracy of major depressive disorders (MDDs), 
identify factors that characterize patients and promote individualized therapy. The goal of this study was to assess a blood-based 
biomarker panel, which showed promise in adolescents with MDD, in adult primary care patients with MDD and age-, gender- and 
race-matched nondepressed (ND) controls. Patients with MDD received cognitive behavioral therapy (CBT) and clinical assessment 
using self-reported depression with the Patient Health Questionnaire-9 (PHQ-9). The measures, including blood RNA collection, 
were obtained before and after 18 weeks of CBT. Blood transcript levels of nine markers of yADC\^3, DGKA, FAM46A, IGSF4A/CADM1 , 
KIAA1539, MARCKS, PSMEl, RAPHl and TLR7, differed significantly between participants with MDD (N = 32) and ND controls (N = 32) 
at baseline {q < 0.05). Abundance of the DGKA, KIAA1539 and RAPHl transcripts remained significantly different between subjects 
with MDD and ND controls even after post-CBT remission (defined as PHQ-9 < 5). The ROC area under the curve for these 
transcripts demonstrated high discriminative ability between MDD and ND participants, regardless of their current clinical status. 
Before CBT, significant co-expression network of specific transcripts existed in MDD subjects who subsequently remitted in 
response to CBT, but not in those who remained depressed. Thus, blood levels of different transcript panels may identify the 
depressed from the nondepressed among primary care patients, during a depressive episode or in remission, or follow and predict 
response to CBT in depressed individuals. 
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INTRODUCTION 

Major depressive disorder (MDD) is a complex psychiatric disease 
affecting 6.7% of the US adult population in a year,^ with a steadily 
increasing prevalence.^'^ The World Health Organization has 
determined that depressive disorders are the leading cause of 
disability worldwide.^ The recurrent and chronic forms of 
depression account for the bulk of the high burden associated 
with the disorder.^ Approximately 50% of patients who recover 
from an initial episode of depression will have at least one 
subsequent depressive episode, and those patients with a history 
of two or more past episodes will have a 70-80% likelihood of 
recurrence in their lives. The number of previous episodes and 
subclinical residual symptoms are the most important predictors 
of relapse.^ Successful treatment of depression relies first upon 
accurate diagnosis. Currently, depression is diagnosed based on 
the patients' self-report of their symptoms and the evaluation of a 
structured psychiatric interview(s) with the patient by a psychia- 
trist, psychologist or primary care physician. The severity of self- 
reported and clinician rated depression symptoms is sometimes in 
disagreement, since depressed patients frequently underreport 
depression symptoms or inadequately characterize them.^ Simi- 
larly, there is no perfect harmony between the different 



depression scoring systems. There are no reliable, objective 
biological diagnostic tests for major depression either.^ The 
diagnostic delay is 2-40 months; and the longer this period, 
more difficult it is to treat depression.^ Most depression is 
treated in primary care: it is estimated that 12.5% of primary care 
patients have had MDD in any given year, but only 47% of those 
cases are recognized clinically.^^ 

Biological markers of different etiologic pathways and/or 
endophenotypes could provide objective, laboratory-based data 
to augment patient reported symptoms and improve the accuracy 
of diagnosis. Biomarkers may be useful in providing information 
about prognosis under different treatment regimens. The 
diagnosis and classification of MDD could lead to enhanced 
individualized care. We have recently identified a panel of blood 
biomarkers in specific animal models of depression^ ^ and showed 
that they can distinguish subjects with early-onset MDD from 
matched controls.^"^ The present study explores these candidate 
biomarkers in a clinical setting, in which subjects with MDD are 
treated with cognitive behavioral therapy (CBT). 
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SUBJECTS AND METHODS 

This study was conducted in the context of a previously reported trial 
comparing the efficacy of face-to-face and telephone-administered CBT.^^ 
Participants with MDD were randomized to two groups, comparing 18 
sessions of CBT, administered either face-to-face or by telephone by a PhD 
level psychologist. These methods of treatment administration produced 
equivalent outcomes for depression.^^ 

Recruitment and procedures 

Participants, recruited from November 2007 to December 2010 were from 
Northwestern General Internal Medicine clinics and were included if they 
were ^ 18 years old, could speak and read English and provided informed 
consent. MDD subjects met criteria for MDD based on the Mini 
International Neuropsychiatry Interview,^ ^ had a Hamilton Depression 
Rating Scale (Ham-D)^^ score ^16 and were able to participate in face-to- 
face or telephone therapy. Nondepressed controls matched by age, race 
and sex, were recruited with the help of NUgene (Northwestern Medicine, 
Northwestern University) and included if they did not meet criteria for 
depression and scored:^ four on Patient Health Questionnaire-9 (PHQ-9). 

Participants were excluded if they: (1) had visual or hearing impairments; 
(2) met diagnostic criteria for other psychiatric disorders (for example, 
bipolar disorder, psychotic disorders, and so on) or depression of organic 
etiology (for example, hypothyroidism); (3) reported alcohol or substance 
abuse severe enough that two psychologists agreed psychotherapy would 
be inappropriate, (4) met criteria for dementia by scoring < 25 on the 
telephone interview for cognitive status;^ ^ (5) exhibited severe suicidality, 
including a plan and intent, or a suicide attempt in the past 5 years; (6) 
were receiving or planning to receive individual psychotherapy; or (7) had 
initiated antidepressant pharmacotherapy in the previous 10 days. This 
trial was approved by the Northwestern University IRB and was monitored 
by an independent data safety monitoring board. 

MDD diagnosis was evaluated at baseline and at posttreatment (week 
18) using the Mini International Neuropsychiatry Interview. Depression 
severity was evaluated by self-report using the PHQ-9 at baseline and 
posttreatment. The PHQ-9 was selected because its items map onto 
diagnostic criteria for MDD and are, therefore, specific to depression. 
Because patients who no longer reach criteria for MDD at posttreatment 
may still experience significant symptoms, patients were also classified as 
reaching full remission based on the MacArthur criterion of PHQ-9 5.^° 
Venous blood (2.5 ml) was collected into PAXgene Blood RNA tubes 
(Qiagen, Germantown, MD, USA) for RNA analysis both at baseline and 
posttreatment from subjects with MDD and at baseline from ND controls. 

Determination of blood transcript levels of candidate markers 
Blood RNA was extracted using the PAXgene Blood RNA Kit (Qiagen), 
according to the manufacturer's protocol. The yield and quality of 
extracted RNA were assessed using the NanoDrop 1000 spectrophot- 
ometer (NanoDrop Technologies, Wilmington, DE, USA). cDNA was 
prepared using random primers and the TaqMan RT reagents (ABI). 
Specific primers are listed in Supplementary Table SI, some of which were 
reported previously.^^ Quantitative real time PGR (qPCR) was carried out 
using SYBR Green and the ABI 7900 (Applied Biosystems, Foster City, CA, 
USA), with 18s rRNA as the internal control. 18s was chosen as the 
normalizing gene after an extensive comparative pilot study using several 
housekeeping genes (data not shown). Each sample for each transcript was 
measured in triplicate using 5 ng of template. 

Baseline samples, collected from 32 patients with a current MDD 
episode and 32 ND subjects, were analyzed in parallel. ACT values from 
qPCR characterized transcript abundance, where ACT is the cycle threshold 
difference between the target gene and the housekeeping gene. Then, to 
avoid inter-assay variations, we re-assayed the triads: baseline MDD, post- 
CBT MDD and ND samples. In the statistical analyses of these results, we 
either compared MDD and ND transcript levels using the ACT values, or 
took advantage of the control matching design and normalized the ACT 
values of MDD subject to that of the matching ND control (which is 1) and 
this relative quantification number (RQ) was used. 

Statistical analyses 

Baseline PHQ-9 and ACT transcript levels were compared between MDD 
and ND using paired f-tests, and presented with mean and standard 
deviation. To avoid inflated type I errors, false discovery rates^^ were 
estimated by q-values using the q-value package^^ in R.^^ We made an a 



priori decision to interpret significant (^-values and an effect size of >0.45 
as clinically significant, which corresponds to Cohen's medium effect size.^^ 

Posttreatment RQ were compared between those who were in clinical 
remission (PHQ-9 :^5) and those who remained depressed using general- 
ized linear models, adjusting for age, sex and race. Nonparametric tests 
were also used with similar results (not reported). Note that when 
comparisons were made between matched data (MDD and ND), analyses 
were not further adjusted for age, sex and race; however, in the 
comparisons of remission or not at week 18, analyses were adjusted. 
Correlational expression analyses of RQ values were carried out by 
Spearman correlation. Two genes were considered to be co-expressed if 
their Spearman correlation was >0.7 and the test statistics against zero 
correlation was <0.01. 

Area under the receiver operating curve (AUG) was calculated against 
ND for those who remitted and those who did not. As an alternative 
analysis,^^ 1000 bootstrap samples of size 50 (including non-remitters, 
remitters and ND) were taken, and AUG was calculated for non-remitters 
relative to ND, and remitters relative to ND, separately. Median AUG and 
empirical 95% intervals are presented for each comparison. Transcripts 
showing significant AUG using both methods of analyses are reported. 

RESULTS 

Participant characteristics 

A total of 64 participants were enrolled, 32 with MDD and 32 who 
were not depressed (ND). As expected, the MDD sample was 
significantly more depressed on the PHQ-9 measure at baseline 
(P< 0.001; Table 1). 

Posttreatment depression scores were available for 27 MDD 
subjects. There was a significant (P< 0.001) reduction in PHQ-9, 
dropping from baseline 18.4 (s.d. = 4.5) to 8.2 (s.d. = 6.1). Of these 
27 MDD participants, 22 completed week 18 blood draws. MDD 
subjects were characterized posttreatment, using the MacArthur 
criteria of PHQ-9 < 5 for full remission,^° into those in full 
remission {N = 9) and those who were not (A/ =13). Patients 
meeting the remission criterion had a mean PHQ-9 = 2.44 (s. 
d.= 1.94) whereas those who did not meet the criterion had a 
mean PHQ-9 = 11.46 (s.d. = 4.71). Patients who remitted were 
younger (age 42.8 (s.d. = 1 1.4) vs 57.4 (s.d.= 14.7), Student's f-test 
P = 0.008), and more likely employed (82% vs 38%, Fisher's exact 
test, P = 0.047). There were no differences between remitters and 
non-remitters in comorbidities, education, marital status, race, sex 
or antidepressant use. 

Transcript level differences at baseline and post-CBT 
Levels of twenty transcripts were measured at both baseline and 
post-CBT. (Tables 2, 3 and Supplementary Figures 1 and 2). At 
baseline, blood transcript levels differed between subjects with 
MDD from their age, sex and race-matched ND controls for 
adenylate cyclase 3 {ADCY3), diacylglycerol kinase, alpha (DGKA), 
family with sequence similarity 46, member A (FAM46A), 
immunoglobulin superfamily, member 4 {IGSF4) also known as 
cell adhesion molecule 1 (CADMl), KIAA1539, which either have no 
brain-related functions to date or no known function, myristoy- 
lated alanine-rich protein kinase C substrate (MARCKS), protea- 
some activator subunit 1 (PSMEl), Ras association and pleckstrin 
homology domains 1 also known as LPD (RAPHl) and intracellular 
Toll-like receptor 7 (TLR7) (q-values < 0.05; Table 2). Post-CBT 
transcript abundance in MDD samples remained significantly 
different from ND controls for DGKA, IGSF4A/CADMI KIAA1539, 
MARCK5 and RAPHl (Table 3). There was no difference in baseline 
gene expression between the group assigned to face-to-face CBT 
compared with telephone-administered CBT (mean Wilcoxon 
rank-sum P = 0.679 (range = 0.071 -0.969)); nor were there differ- 
ences observed by age, sex, education, marital status, race, 
antidepressant medication use, employment, history of general- 
ized anxiety, history of obsessive compulsive disorder, alcohol 
dependence, melancholia, panic disorder or social phobia. 
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Table 2. Baseline transcript levels (ACT) for subjects with MDD 
(A/ = 32) and sex-, age-, race-matched ND controls (A/ =32) 



Gene 


MDD mean 


CONTROL 


P-value 


q 


Cohen's 




fs.d.) 


mean (s.d.) 








ADCY3 


20.95 (2.07) 


21.31 (2.18) 


0.012 


0.026 


0.473 


AMFR 


19.28 (0.87) 


19.27 (0.42) 


0.995 


0.995 


0.001 


ASAH1 


16.74 (2.45) 


16.88 (2.25) 


0.352 


0.503 


0.167 


ATPnC 


19.22 (1.83) 


19.20 (2.17) 


0.902 


0.949 


0.022 


CD59 


19.98 (0.65) 


19.89 (0.66) 


0.550 


0.651 


0.107 


CDR2 


20.60 (1.13) 


20.49 (0.65) 


0.554 


0.652 


0.106 


CMAS 


22.41 (1.36) 


22.98 (0.69) 


0.046 


0.084 


0.367 


DGKA 


16.69 (0.80) 


17.41 (0.75) 


0.000 


0.003 


0.730 


FAM46A 


19.65 (2.27) 


19.96 (2.40) 


0.005 


0.017 


0.532 


IGSF4A/ 


22.22 (3.07) 


23.05 (2.57) 


0.002 


0.010 


0.597 


CADM1 












KIAA1539 


18.22 (0.61) 


18.64 (0.45) 


0.001 


0.009 


0.625 


MAF 


22.00 (2.31) 


22.08 (2.47) 


0.637 


0.708 


0.084 


MARCKS 


16.85 (2.31) 


17.41 (2.22) 


0.005 


0.017 


0.533 


NAGA 


1 7.63 (2.22) 


1 7.85 (2.30) 


0.041 


0.083 


0.376 


PSME1 


17.17 (0.75) 


1 7.55 (0.52) 


0.008 


0.023 


0.501 


PTP4A3 


20.79 (0.82) 


21.03 (0.62) 


0.154 


0.257 


0.258 


RAPH1 


21.73 (1.47) 


22.80 (0.54) 


< 0.001 


0.003 


0.737 


SLC4A1 


1 7.08 (0.85) 


1 7.22 (0.85) 


0.454 


0.606 


0.134 


TLR7 


19.53 (2.36) 


19.89 (2.35) 


0.009 


0.023 


0.493 


ZNF291 


23.55 (3.18) 


23.20 (3.20) 


0.176 


0.270 


0.253 



Abbreviations: MDD, major depressive disorder; ND, nondepressed. Bold 
entries represent significant results. 



Table 3. 



(A/ = 22) and sex-, age-, race-matched ND controls (A/ =22) 


Gene 


MDD mean 


CONTROL 


P-value 


q 


Cohen's D 




fs.d.) 


mean (s.d.) 








ADCY3 


20.25 (1.85) 


20.44 (1.96) 


0.092 


0.190 


0.357 


AMFR 


19.33 (0.75) 


19.33 (0.47) 


0.913 


0.913 


0.004 


ASAHl 


15.98 (2.15) 


16.04 (2.13) 


0.680 


0.716 


0.065 


ATPllC 


18.78 (2.13) 


18.48 (2.11) 


0.123 


0.208 


0.363 


CD59 


20.09 (0.67) 


19.79 (0.74) 


0.125 


0.208 


0.355 


CDR2 


20.68 (1.23) 


20.46 (0.66) 


0.376 


0.470 


0.215 


CMAS 


22.55 (1.04) 


23.04 (0.70) 


0.076 


0.190 


0.481 


DGKA 


16.33 (0.72) 


1 7.40 (0.89) 


< 0.001 


0.005 


0.883 


FAM46A 


18.84 (2.18) 


19.04 (2.22) 


0.232 


0.331 


0.309 


IGSF4A/ 


21.50 (2.48) 


22.19 (2.48) 


0.001 


0.005 


0.862 


CADM1 












KIAA1539 


18.15 (0.62) 


18.71 (0.50) 


< 0.001 


0.005 


0.840 


MAF 


21.16 (2.11) 


21.10 (2.33) 


0.555 


0.653 


0.026 


MARCKS 


16.22 (1.82) 


16.64 (2.12) 


0.010 


0.040 


0.558 


NAGA 


16.89 (2.15) 


16.99 (2.14) 


0.595 


0.661 


0.119 


PSMEl 


1 7.04 (0.68) 


1 7.46 (0.56) 


0.033 


0.094 


0.533 


PTP4A3 


20.79 (0.74) 


21.05 (0.64) 


0.163 


0.251 


0.339 


RAPH1 


21.45 (1.07) 


22.86 (0.60) 


< 0.001 


0.005 


1.086 


SLC4A1 


17.02 (1.16) 


1 7.37 (0.80) 


0.095 


0.190 


0.290 


TLR7 


18.52 (2.24) 


18.92 (2.09) 


0.016 


0.533 


0.550 


ZNF291 


22.48 (3.12) 


21.95 (2.96) 


0.373 


0.470 


0.249 



Abbreviations: MDD, major depressive disorder; ND, nondepressed. Bold 
entries represent significant results. 



Scatter plots of ACT values for each transcript measured at 
baseline and post-CBT in MDD subjects and their nnatched 
controls are presented in Supplementary Figures 1 and 2. 

Transcript level differences by remission 

Subjects were divided into those who remitted at the end of CBT 
(PHQ-9<5), and those who did not remit. Post-CBT RQ values 
differed between these subjects in levels of N-acylsphingosine 
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Table 4. Posttreatment transcript levels (RQ; normalized to matched 
ND controls) of MDD patients by remission status 



Gene 


PHQ-9>5 


PHQ-9<5 


Age-, sex-, 


Co lien's D 




(M = 13) RQ 


(M = 9) RQ 


race- 






mean (5E) non- 


mean (SE) 


adjusted 






remitted 


remitted 


P-value 




ADCY3 


1.52 (0.22) 


1.47 (0.22) 


0.879 


0.033 


AMFR 


1.33 (0.24) 


1.47 (0.25) 


0.669 


0.093 


ASAH1 


0.90 (0.12) 


1.38 (0.12) 


0.009 


0.630 


ATPnC 


1.12 (0.16) 


0.62 (0.16) 


0.035 


0.486 


CD59 


0.92 (0.15) 


1.13 (0.16) 


0.316 


0.221 


CDR2 


1.54 (0.26) 


1.22 (0.27) 


0.377 


0.193 


CMAS 


2.13 (0.64) 


2.69 (0.68) 


0.523 


0.139 


DGKA 


3.56 (0.92) 


3.78 (0.97) 


0.863 


0.037 


FAM46A 


1.41 (0.23) 


1.20 (0.23) 


0.518 


0.141 


IGSF4A/ 


1.60 (0.44) 


2.08 (0.44) 


0.433 


0.171 


CADMl 










KIAA1539 


1.72 (0.21) 


2.34 (0.22) 


0.042 


0.469 


MAF 


1.08 (0.21) 


1.09 (0.21) 


0.963 


0.010 


MARCKS 


1.26 (0.33) 


2.18 (0.33) 


0.058 


0.432 


NAGA 


1.25 (0.21) 


1.28 (0.21) 


0.910 


0.023 


PSMEl 


1.70 (0.30) 


2.30 (0.31) 


0.148 


0.323 


PTP4A3 


1.52 (0.25) 


1.95 (0.26) 


0.206 


0.280 


RAPHl 


4.30 (1.29) 


5.49 (1.36) 


0.504 


0.146 


SLC4A1 


1.70 (0.52) 


2.47 (0.54) 


0.279 


0.238 


TLR7 


1.77 (0.29) 


1.76 (0.29) 


0.979 


0.006 


ZNF291 


1.39 (0.44) 


1.21 (0.44) 


0.771 


0.063 



Abbreviations: MDD, major depressive disorder; ND, nondepressed; RQ, 
relative quantification number. Bold entries represent significant results. 



annldohydroiase (acid cerannidase) 1 (ASAHl), ATPase, Class VI, 
Type lie (ATPllC) and KIAA1539 (Table 4). These differences, 
however, were not significant when adjusted for nnultiple 
connparisons (q-values not shown). Post-CBT, there was no 
difference in transcript levels between the group assigned to 
face-to-face CBT connpared with telephone-adnninistered CBT 
(nnean Wilcoxon rank-sunn P= 0.694 (range = 0.336-0.999)). 

Trait nnarkers: transcript level differences between subjects with 
MDD and controls regardless of current clinical diagnosis 
Both baseline and post-CBT blood transcript levels differed 
between MDD vs ND controls in DGKA, IGSF4A/CADI\/n , KIAA1539, 
MARCKS and RAPHl. Rec eiver-operating characteristics for post- 
CBT nneasures of blood transcript levels revealed that RAPHl, 
KIAA1539 and DGKA had significant AUG for subjects who rennitted 
after CBT, but also for subjects who did not remit (Supplennentary 
Table S2). Thus, these nnarkers had good discrinninative character- 
istics, with significant specificity and selectivity for subjects with 
MDD vs ND controls even when MDD patients were no longer 
depressed. 

Co-expression patterns of transcripts at baseline predict rennission 
or persisting clinical depression post CBT 

Correlational expression networks representing baseline co- 
expression patterns of transcripts, normalized to their matched 
controls (RQ), in subjects remitted and non-remitted after CBT are 
shown in Figure 1. In total, 15 co-expressed gene-pairs were found 
in the non-remitted group and 20 pairs were found in the 
remitted group. The extensiveness of correlation at baseline was 
greater for those who remitted compared with those who did not 
at post-CBT. Eleven 'hubs', identified as having significant 
correlations with three or more transcripts, were in the remitted, 
compared with five in the non-remitted population of MDD 
patients. Specifically, transcript levels of ADCY3, DGKA, IGSF4A/ 
CADMl, PSMEl and RAPHl are correlated significantly at baseline 
in subjects who remitted in response to CBT. None of these 
transcripts correlated in their expression at baseline in those 
subjects who remained clinically depressed after CBT. 



DISCUSSION 

The search for an objective diagnostic test with high sensitivity 
and specificity for MDD has been ongoing for decades. We and 
others have shown that it is possible to approach this question 
using unbiased 'omics' to identify candidate blood marker panels, 
as reviewed by Menke.^^ In the present work, we examined the 
performance of a panel of blood transcriptomic markers in a 
clinical population of subjects with MDD. The results are the first at 
many levels. Abundance of nine transcripts differed significantly 
between subjects with MDD and ND controls, suggesting that this 
panel could be diagnostic in a clinical, not medication-free 
population. Transcript levels of three genes follow the remission 
status post CBT, raising the possibility that these can be used to 
survey treatment efficacy. Blood levels of three transcripts differed 
between ND controls and MDD subjects, in those who remitted in 
response to CBT, and in those who did not respond to CBT and 
remained clinically depressed. These markers come closest to the 
ultimate goal of identifying predisposition to depression, even in 
the absence of a current depressive episode. Finally, another first 
is the fingerprint-like co-expression network that could predict 
whether the patient will get well after CBT or not, before the 
therapy begins. 

In our previous studies, we described the discovery of a blood 
transcriptomic marker panel using two animal models of 
depression, a genetic and a chronic stress model. ^^'^^ In a pilot 
study we had shown that some of the transcriptomic markers 
comprising this panel differ in their abundance between 
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teenagers with MDD and those with no disorder.^"^ The present 
study differs in the wide age range (23-83) and concurrent 
medication status of the participants compared with the tight age 
range (15-19) and medication-free requirement of the adolescent 
study. Despite these differences, three of the nine transcripts, 
whose blood levels differed significantly between subjects with 
MDD and ND controls in the current study, PSMEl, RAPHl and 
IGSF4A also differed between medication-free MDD teenagers and 
controls.^"^ These results, therefore, represent a confirmation and 
suggest that these transcripts participate in processes that are 
characteristic of depressed mood, or of other endophenotypes 
intricately involved in MDD. 

Earlier studies of transcriptional profiling of psychiatric illnesses 
used RNA derived from postmortem brain tissue.^'' However, the 
limitation of this approach and accessibility of peripheral blood 
cells prompted a series of investigations to prove the overlap 
between gene expression in peripheral blood and brain.^^"^° 
Transcriptional profiling in peripheral blood has been used to find 
biomarkers for patients with psychiatric disorders.^^"^^ A multi- 
assay, serum-based biological diagnostic tests for MDD has also 
been proposed by Papakostas et al?^ Serum levels of nine 
biomarkers are measured in this test and a proprietary algorithm 
generates the results. These blood-based biomarker studies, 
together with our present and previous^^ reports, demonstrate 
the potential of using blood as a proxy for the brain.^^'^° 

Here we report the first identification of molecular trait markers 
for MDD. Three transcripts, with levels that differed significantly 
between MDD patients and controls when subjects with MDD 
were depressed, remained different even when patients were no 
longer clinically depressed. This pattern may be due to direct or 
indirect consequences of genetic variation between controls and 
MDD subjects. However, none of these three genes, RAPHl, which 
is also named as amyotrophic lateral sclerosis 2 (juvenile) 
chromosome region, candidate 9, KIAA1539, also known as 
FAM214B, and DGKA were found to be associated with major 
depression in individual genome-wide association studies or 
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Figure 1. Correlational expression networks. Co-expression patterns of transcript levels at baseline, nornnalized to their nnatched controls (RQ), 
of subjects rennitted (a) and non-rennitted (b) after CBT Two genes were considered to be co-expressed if their Spearnnan correlation was >0.7 
and the test statistics against zero correlation was < 0.01. Red lines represent connnnon co-expression of non-rennitted and rennitted groups, 
and black lines represent group-specific co-expression. Transcripts showing significantly different levels at baseline between subjects with 
MDD and ND controls are plotted with filled circles. CBT, cognitive behavioral therapy; MDD, nnajor depressive disorder; ND, nondepressed. 



nneta-analyses.^^"'^° Thus, sequence variation or copy number 
variation in other functionally connected genes nnight be 
responsible for the observed trait-marker status. The simultaneous 
transcript level differences from ND controls in these three genes 
could possibly mark subjects with future recurrent depression. 
Liotti et alf^^ used a similar design to identify depression-specific 
cerebral blood flow changes. They found that the pattern of acute 
flow changes in unipolar euthymic patients in full remission is 
distinct from that seen in euthymic healthy volunteers and mirrors 
the untreated depressed state seen during a major depressive 
episode. Thus, our findings along with others, suggest that 
disease-specific modifications of pathways present in unipolar 
depression are independent of clinical illness status, providing 
implications for the vulnerability of patients for illness relapse. 
Should these vulnerability markers be confirmed in larger studies, 
they could distinguish MDD patients from ND controls, with or 
without current active depressive episode. These patients, who are 
likely to have more than one depressive episode, once identified 
would benefit from active long-term antidepressant or psy- 
chotherapy treatment, which could protect, at least in part, from 
relapse."^^ 

After CBT, blood transcript levels of ASAHh ATPllC and 
KIAA1539 differed between remitted and non-remitted MDD 
patients. In a recent study, four transcripts, different from the 
present ones, were identified as peripheral transcriptional markers 
of responders and non-responders to differing antidepressant 
treatment.^^ These markers, therefore, may be markers of 
depressive state, as opposed to an underlying vulnerability. Our 
confidence in the relationship between the presently identified 
markers described here and depressive state is strengthened by 
the fact that in our study the treatment was psychological. 
Furthermore, as most remitted and non-remitted patients were 
medicated (Supplementary Table S3), the medications themselves 
did not interfere with the blood marker measurements. 

In a comprehensive and timely summary, Gaiteri et alf^ suggest 
that co-expression networks may reflect endogenous regulatory 



systems, and co-expression modules can be significant predictors 
of disease. Moreover, they propose that disease severity may 
interact with molecular connectivity in such a way that it places 
disease genes at different network locations. This possibility is 
explored here when we found that differences in co-expression 
patterns could become a predictor of CBT efficacy before therapy. 
Specifically, concomitant differences in the abundance of ADCY3, 
DGKA, IGSF4A/ CADMh PSMEl and RAPH1 transcripts between 
MDD subjects and ND controls could have predictive value. This is 
a very promising finding as to date, there are no biomarkers or 
predictors of effectiveness of psychotherapy. Biomarkers of 
antidepressant therapy employed pharmacogenomics studies 
with modest results.^^'"^^ A nonmolecular frontal quantitative 
electroencephalographic biomarker has also been proposed as a 
predictor of response to escitalopram after 1 week of treatment."^^ 
Thus, the co-expression criteria of the specific markers before 
treatment could be developed, if confirmed, into an aid of CBT 
treatment choice in clinical practice. A current meta-analysis 
identified that only 14% more patients achieve remission after 
psychotherapy compared with care-as-usual,^^ but this number 
may be substantially greater if patients who are likely to respond 
to CBT could be identified ahead of time. The significance of this 
possibility is clearly stated by Licinio and Wong"^^'"^^ who state that 
when we discover a priori biomarkers, identifying which patients 
will respond to which drugs (or treatment, as we suggest), 
therapeutics in psychiatry will be much further ahead than it 
is now. 

This study has a number of limitations. These include the 
modest sample size. In addition, there were only two samples for 
each MDD subject, and only one ND sample. Repeated sampling 
of the ND subjects would have added to the power of the design, 
although gene expression patterns are thought to be stable 
within, but differ between individuals.^^"" Should we have 
measured levels of transcripts at different time points during 
and after CBT, the kinetics of changes in transcript levels could 
have been determined in parallel with those of mood. An 
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Interesting finding that points to the need of large number of 
age-, gender- and race-stratified controls was that baseline 
gene expression of MDD subjects and their matched controls 
did vary between pairs with and without a second MDD draw 
(Supplementary Table S4). The cause of this is unknown, although 
there seems to be an age difference between women who did and 
did not give blood at post-CBT (see Table 1) and, therefore, an age 
difference in the matching control sample. This age difference did 
not reach significance, due to the very large age range of our 
sample, but it could affect gene expression differences. Even more 
interesting is that although matching was done only by age, sex 
and race, gene expression overall was not strongly associated with 
any of these factors. Still, transcript levels in the blood of both 
MDD and ND subjects varied at baseline by status of the second 
draw. It is worthwhile to note that the 'trait' marker candidates, 
DGKA, KIAA1359 and RAPHl showed none of the differences 
between the groups with and without post-CBT blood sample. 

CONCLUSIONS 

This study identified nine transcripts that differentiated depressed 
from nondepressed participants, three of which can distinguish 
subjects with MDD from nondepressed controls, even after 
remission. A three-transcript panel discerned remitted from non- 
remitted patients post CBT, and concomitant abundance differ- 
ences in some transcripts at baseline might indicate treatment 
responsiveness a priori. Thus, blood levels of different transcript 
panels may be useful in identifying depressed primary care 
patients, as well as in predicting and monitoring response to CBT. 
Future studies are aimed at validating these results in a larger 
patient population and determining the selectivity of these panels 
in patient population that include other psychiatric illnesses. 

CONFLICT OF INTEREST 

EER is named as an inventor on two pending patent applications, filed and owned by 
Northwestern University. The remaining authors declare no conflict of interest. 

ACKNOWLEDGMENTS 

This study was supported by the National Institute of Mental Health (grant R21 
MH077234 and R01 MH059708) and by grants from the Davee Foundation. The 
sponsors had no role in the design and conduct of the study; collection, 
management, analysis or interpretation of the data; preparation, review or approval 
of the manuscript; and decision to submit the manuscript for publication. We 
acknowledge Neha Mehta for her help in finalizing the manuscript. We thank all 
subjects for the participation in this research. 

REFERENCES 

1 Kessler RC, Berglund P, Demler O, Jin R, Merikangas KR, Walters EE. Lifetime 
prevalence and age-of-onset distributions of DSM-IV disorders in the National 
Comorbidity Survey Replication. Arch Gen Psychiatry 2005; 62: 593-602. 

2 Hidaka BH. Depression as a disease of modernity: explanations for increasing 
prevalence. J Affect Disord 2012; 140: 205-214. 

3 Lepine JP, Briley M. The increasing burden of depression. Neuropsychiatr Dis Treat 
2011; 7: 3-7. 

4 Ferrari AJ, Charlson FJ, Norman RE, Patten SB, Freedman G, Murray CJ et al. Burden 
of depressive disorders by country, sex, age, and year: findings from the global 
burden of disease study 2010. PLoS Med 2013; 10: el 001 547. 

5 Hardeveld F, Spijker J, De Graaf R, Nolen WA, Beekman AT. Prevalence and pre- 
dictors of recurrence of major depressive disorder in the adult population. Acta 
Psychiatr Scand 2010; 122: 184-191. 

6 ten Doesschate MC, Bockting CL, Koeter MW, Schene AH. Prediction of recurrence 
in recurrent depression: a 5.5-year prospective study. J Clin Psychiatry 2010; 71: 
984-991. 

7 Carter JD, Frampton CM, Mulder RT, Luty SE, Joyce PR. The relationship of 
demographic, clinical, cognitive and personality variables to the discrepancy 
between self and clinician rated depression. J Affect Disord 2010; 124: 202-206. 



8 Riedel M, Moller HJ, Obermeier M, Schennach-Wolff R, Bauer M, Adii M et al. 
Response and remission criteria in major depression~a validation of current 
practice. J Psychiatr Res 2010; 44: 1063-1068. 

9 Belmaker RH, Agam G. Major depressive disorder. N Engl J Med 2008; 358: 55-68. 

10 Altamura AC, Buoli M, Albano A, Dell'Osso B. Age at onset and latency to treat- 
ment (duration of untreated illness) in patients with mood and anxiety disorders: 
a naturalistic study. Int Clin Psychopharmacol 2010; 25: 172-179. 

11 Huerta-Ramirez R, Bertsch J, Cabello M, Roca M, Haro JM, Ayuso-Mateos JL. 
Diagnosis delay in first episodes of major depression: a study of primary care 
patients in Spain. 7 /Affecf Disord 2013; 150: 1247-1250. 

1 2 Pence BW, O'Donnell JK, Gaynes BN. The depression treatment cascade in primary 
care: a public health perspective. Curr Psychiatry Rep 2012; 14: 328-335. 

13 Andrus BM, Blizinsky K, Vedell PT, Dennis K, Shukia PK, Schaffer DJ et al. Gene 
expression patterns in the hippocampus and amygdala of endogenous depres- 
sion and chronic stress models. Mol Psychiatry 2012; 17: 49-61. 

14 Pajer K, Andrus BM, Gardner W, Lourie A, Strange B, Campo J et al. Discovery of 
blood transcriptomic markers for depression in animal models and pilot valida- 
tion in subjects with early-onset major depression. TransI Psychiatry 2012; 2: el 01. 

15 Mohr DC, Ho J, Duffecy J, Reifler D, Sokol L, Burns MN et al. Effect of telephone- 
administered vs face-to-face cognitive behavioral therapy on adherence to 
therapy and depression outcomes among primary care patients: a 
randomized trial. JAMA 2012; 307: 2278-2285. 

16 Mohr DC, Vella L, Hart SL, Heckman T, Simon GE. The effect of telephone- 
administered psychotherapy on depression and attrition: a meta-analysis. Clin 
Psych 2008; 1 5: 243-253. 

17 Sheehan DV, Lecrubier Y, Sheehan KH, Amorim P, Janavs J, Weiller E et al. The 
Mini-International Neuropsychiatric Interview (M.I.N. I.): the development and 
validation of a structured diagnostic psychiatric interview for DSM-IV and ICD-10. 
J Clin Psychiatry 1998, 59(Suppl 20): 22-33. quiz 34-57. 

18 Hamilton M. A rating scale for depression. J Neurol Neurosurg Psychiatry 1960; 23: 
56-62. 

19 Desmond DW, Tatemichi TK, Hanzawa L. The Telephone Interview for Cognitive 
Status (TICS): reliability and validity in a stroke sample. IntJGeriatr Psychiatry ^994•, 
9: 803-807. 

20 MacArthur Foundation's Initiative on Depression and Primary Care. The MacArthur 
Initiative on Depression and Primary Care at Dartmouth and Duke: Depression 
Management Toolkit. Hanover, NH, USA: Dartmouth, 2004. 

21 Benjamini Y, Hochberg Y. Controlling the False Discovery Rate - a Practical and 
Powerful Approach to Multiple Testing. J Roy Stat Soc B Met 1995; 57: 289-300. 

22 Storey JD, Tibshirani R. Statistical significance for genomewide studies. Proc Natl 
Acad Sci USA 2003; 100: 9440-9445. 

23 Team RC R: A language and environment for statistical computing. ISBN 
3-900051-07-0. R Foundation for Statistical Computing. Vienna, Austria, 2013. 
Available online at http://www.R-project.org; 2005. 

24 Cohen J. Statistical Power Analysis for the Behavioral Sciences. 2nd edn, Hillsdale, 
NJ, USA: L Eribaum Associates, 1988. 

25 Pepe M, Longton G, Janes H. Estimation and Comparison of Receiver Operating 
Characteristic Curves. Stata J 2009; 9: 1. 

26 Menke A. Gene expression: biomarker of antidepressant therapy? Int Rev Psy- 
chiatry 201 3; 25: 579-591 . 

27 Mehta D, Menke A, Binder EB. Gene expression studies in major depression. Curr 
Psychiatry Rep 2010; 12: 135-144. 

28 Liew CC, Ma J, Tang HC, Zheng R, Dempsey AA. The peripheral blood tran- 
scriptome dynamically reflects system wide biology: a potential diagnostic tool. J 
Lab Clin Med 2006; 147: 126-132. 

29 Rollins B, Martin MV, Morgan L, Vawter MP. Analysis of whole genome biomarker 
expression in blood and brain. Am J Med Genet B Neuropsychiatr Genet 2010; 
153B: 919-936. 

30 Sullivan PF, Fan C, Perou CM. Evaluating the comparability of gene expression in 
blood and brain. Am J Med Genet B Neuropsychiatr Genet 2006; 141 B: 261-268. 

31 Kurian SM, Le-Niculescu H, Patel SD, Bertram D, Davis J, Dike C et al. Identification 
of blood biomarkers for psychosis using convergent functional genomics. Mol 
Psychiatry 2011; 16: 37-58. 

32 Le-Niculescu H, Kurian SM, Yehyawi N, Dike C, Patel SD, Edenberg HJ et al. 
Identifying blood biomarkers for mood disorders using convergent functional 
genomics. Mol Psychiatry 2009; 14: 156-174. 

33 Padmos RC, Hillegers MH, Knijff EM, Vonk R, Bouvy A, Staal FJ et al. A dis- 
criminating messenger RNA signature for bipolar disorder formed by an aberrant 
expression of inflammatory genes in monocytes. Arch Gen Psychiatry 2008; 65: 
395-407. 

34 Menke A, Arloth J, Putz B, Weber P, Klengel T, Mehta D et al. Dexamethasone 
stimulated gene expression in peripheral blood is a sensitive marker for gluco- 
corticoid receptor resistance in depressed patients. Neuropsychopharmacology 
2012; 37: 1455-1464. 



Translational Psychiatry (2014), 1 -7 



© 2014 Macmillan Publishers Limited 



Blood transcript markers of major depression 
EE Redei et al 



35 Segman RH, Goltser-Dubner T, Weiner I, Canetti L, Galili-Weisstub E, Milwidsky A 
et al. Blood mononuclear cell gene expression signature of postpartum depres- 
sion. Mol Psychiatry 2010; 15: 2. 

36 Spijker S, Van Zanten JS, De Jong S, Penninx BW, van Dyck R, Zitman FG et al. 
Stimulated gene expression profiles as a blood marker of major depressive dis- 
order. Biol Psychiatry 2010; 68: 179-186. 

37 Papakostas Gl, Shelton RC, Kinrys G, Henry ME, Bakow BR, Lipkin SH et al. 
Assessment of a multi-assay, serum-based biological diagnostic test for major 
depressive disorder: a pilot and replication study. Mol Psychiatry 2013; 18: 
332-339. 

38 Shyn SI, Shi J, Kraft JB, Potash JB, Knowles JA, Weissman MM et al. Novel loci for 
major depression identified by genome-wide association study of Sequenced 
Treatment Alternatives to Relieve Depression and meta-analysis of three studies. 
Mol Psychiatry 2011; 16: 202-215. 

39 Cross-Disorder Group of the Psychiatric Genomics Consortium. Identification of 
risk loci with shared effects on five major psychiatric disorders: a genome-wide 
analysis. Lancet 2013; 381: 1371-1379. 

40 Wray NR, Pergadia ML, Blackwood DHR, Penninx BWJH, Gordon SD, Nyholt DR 
et al. Genome-wide association study of major depressive disorder: new results, 
meta-analysis, and lessons learned. Mol Psychiatry 2012; 17: 36-48. 

41 Liotti M, Mayberg HS, McGinnis S, Brannan SL, Jerabek P. Unmasking disease- 
specific cerebral blood flow abnormalities: mood challenge in patients with 
remitted unipolar depression. Am J Psychiatry 2002; 159: 1830-1840. 

42 Segal ZV, Bieling P, Young T, MacQueen G, Cooke R, Martin L et al. Antidepressant 
monotherapy vs sequential pharmacotherapy and mindfulness-based cognitive 
therapy, or placebo, for relapse prophylaxis in recurrent depression. Arch Gen 
Psychiatry 2010; 67: 1256-1264. 

43 Belzeaux R, Bergon A, Jeanjean V, Loriod B, Formisano-Treziny C, Verrier L et al. 
Responder and nonresponder patients exhibit different peripheral transcriptional 
signatures during major depressive episode. TransI Psychiatry 2012; 2: el 85. 

44 Gaiteri C, Ding Y, French B, Tseng GC, Sibille E. Beyond modules and hubs: the 
potential of gene coexpression networks for investigating molecular mechanisms 
of complex brain disorders. Genes Brain Behav 2014; 13: 13-24. 



45 Gunther EC, Stone DJ, Gerwien RW, Bento P, Heyes MP. Prediction of clinical drug 
efficacy by classification of drug-induced genomic expression profiles in vitro. 
Proc Natl Acad Sci USA 2003; 100: 9608-9613. 

46 Leuchter AF, Cook lA, Marangell LB, Gilmer WS, Burgoyne KS, Howland RH et al. 
Comparative effectiveness of biomarkers and clinical indicators for predicting 
outcomes of SSRI treatment in major depressive disorder: results of the BRITE- 
MD study. Psychiatry Res 2009; 169: 124-131. 

47 Cuijpers P, Karyotaki E, Weitz E, Andersson G, Hollon SD, van Straten A. The effects 
of psychotherapies for major depression in adults on remission, recovery and 
improvement: a meta-analysis. J Affect Disord 2014, 159: 118-126. 

48 Licinio J, Wong ML. Launching the 'war on mental illness'. Mol Psychiatry 2014; 19: 
1-5. 

49 Licinio J, Wong ML. Pharmacogenomics of antidepressant treatment effects. 
Dialogues Clin Neurosci 2011; 13: 63-71. 

50 Radich JP, Mao M, Stepaniants S, Biery M, Castle J, Ward T et al. Individual-specific 
variation of gene expression in peripheral blood leukocytes. Genomics 2004; 83: 
980-988. 

51 Eady JJ, Wortley GM, Wormstone YM, Hughes JC, Astley SB, Foxall RJ et al. Var- 
iation in gene expression profiles of peripheral blood mononuclear cells from 
healthy volunteers. Physiol Genomics 2005; 22: 402-411. 

52 Mohr S, Liew C-C. The peripheral-blood transcriptome: new insights into disease 
and risk assessment. Trends Mol Med 2007; 13: 422-432. 



©®®© 



This work is licensed under a Creative Commons Attribution- 
NonCommercial-NoDerivs 3.0 Unported License. The images or 
other third party material in this article are included in the article's Creative 
Commons license, unless indicated otherwise in the credit line; if the material is not 
included under the Creative Commons license, users will need to obtain permission 
from the license holder to reproduce the material. To view a copy of this license, visit 
http://creativecommons.Org/licenses/by-nc-nd/3.0/ 



Supplementary Information accompanies the paper on the Translational Psychiatry website (http://www.nature.com/tp) 



© 2014 Macmillan Publishers Limited 



Translational Psychiatry (2014), 1 -7 



